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We review the main results we have obtained in the area of high-energy elastic hadron scattering 
and presented in this series of Workshops on Hadron Interactions. After an introduction to some 
basic experimental and theoretical concepts, we survey the results reached by means of four ap- 
proaches: analytic models, model-independent analyses, eikonal models and nonperturbative QCD. 
Some of the ongoing researches and future perspectives are also outlined. 
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I. INTRODUCTION 

"QCD nowadays has a split personality. It 
embodies hard and soft physics, both being 
hard subjects and the softer the harder." 

Yuri Dokshitzer (2001) [H 

Despite the great success of QCD as the field theory 
of hadronic interactions, there still remains some open 
questions and one of them is related to the hadron-hadron 
scattering at high- energies and small momentum transfer 
(soft diffraction). 

The region of high energies is characterized by scatter- 
ing of particles with center of mass energy ^/s > 10 GeV 
^ 10 rrip (the proton mass). From the experimental point 
of view, diffractive processes are associated with a slow 
increase of the total cross sections, the diffraction pat- 
tern in the differential cross section, and rapidity gaps 
in the plots of pseudo rapidity versus azimuthal angle. 
In the theoretical context, diffraction means that the ini- 
tial and final states in the scattering process have the 
same quantum numbers and, therefore, the exchanged 
"object" has the vacuum quantum numbers (Pomeron). 
The soft diffractive processes are generally classified as 
double diffraction dissociation, single diffraction dissoci- 
ation and elastic scattering. Introductoryreviews on the 
area can be found in Refs. 0, i, S S H Q • 

High-energy elastic hadron scattering is the simplest 
soft diffractive process and, at the same time a topi- 
cal problem in high-energy physics. Being associated 
with long distance phenomena perturbative QCD can 
not be applied. On the other hand, the standard non- 
perturbative approach starts with the ground state (vac- 
uum), proceeds with bound states (mesons, barions) and 
eventually reaches the scattering states. However, it 
is obvious that the vacuum is a non-trivial problem. 
Moreover, even assuming some vacuum concept, to treat 
only one gluon field it is necessary to take into account 
more than 30 invariants, and all that becomes a typical 
problem of statistical physics, with specific technical ap- 
proaches, such as Monte Carlo simulation (lattice QCD). 
Although bound states may be described, the point is 
that, presently, we do not know how to calculate elastic 
scattering amplitudes from a pure nonperturbative QCD 
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formalism. 

At this stage, phenomenology certainly plays an im- 
portant role in the search for connections between exper- 
imental data, model descriptions, and the possible devel- 
opment of new calculational schemes in the underlying 
theory (QCD). Here, however, we are faced with another 
kind of problem, namely, the wide variety of model de- 
scriptions, based on different ideas and approaches, not 
always giving enough support for the development of 
novel calculational schemes well founded on QCD. 

Based on the above facts, our main strategy in the in- 
vestigation of the elastic sector is to search for model 
independent information that may be extracted from the 
experimental data, through approaches that have well es- 
tablished bases on the General Principles, theorems and 
bounds from axiomatic quantum field theory (the ana- 
lytic approach) . Simultaneously, we attempt to construct 
phenomenological models, in agreement with the above 
Principles and connected, in some way, with the under- 
lying dynamics of QCD. 

In this review, it is presented some results we have 
obtained in the area of elastic scattering in the last 
years, with focus on high-energy proton-proton (pp) 
and antiproton-proton (pp) elastic scattering. The 
manuscript is organized as follows. In Sec. II we re- 
call some basic experimental and theoretical concepts, 
defining also our notation. In Sees. Ill, IV, V, and VI 
we present the main results we have obtained through- 
out the analytic approach, model independent analyses, 
eikonal models and nonperturbative QCD, respectively. 
In Sec. VII we discuss some perspectives in the area, 
from both experimental and theoretical points of view. 
A summary and some final remarks are the contents of 
Sec. VIII. 



II. BASIC CONCEPTS 

In this section we recall the physical quantities that 
characterize the elastic scattering and shortly review 
some principles, high-energy theorems, and the main for- 
mulas associated with two basic pictures, usually referred 
as s-channel (geometrical/optical picture) and t-channel 
(exchange picture) 0, |4 ^ J(| • 



A. Physical Quantities 

In elastic scattering, the connection between experi- 
mental data and theory is done by means of the invariant 
scattering amplitude, expressed in terms of two Mandel- 
stam variables, generally the center-of-mass (cm.) en- 
ergy squared s and the four-momentum transfer squared 
t = —q^: F = F{s,t). It is expected that spin effects 
decrease as the energy increases (for some recent results 
see 0), and neglecting spin, the physical quantities that 
characterize the elastic scattering process are the differ- 
ential cross section, 



dc , . TT . , , , n 

-(M) = ^|F(M)P, 



(1) 



where k is the cm. momentum, the elastic integrated 
cross section. 



(Tel[S) = J —{s,t)dt, 
the total cross section (Optical Theorem), 
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crtot{s) = — lmi^(s,0), 
the inelastic cross section 

CTmei(s) = CFtot{s) - (Tel{s), 

the p parameter. 



(2) 



P{s) 



ReF(s,0) 
ImF(s,0)' 



(3) 



and the slope parameter. 



B{s) 



d_ 

dt 



, da , , 



(4) 



The corresponding experimental data have been ana- 
lyzed and compiled by the Particle Data Group and can 
be found in Ref. |9( and quoted references. In what 
follows we shall be mainly interested in pp and pp data 
in the regions: 13.8 GeV < \/s < 1.8 TeV and 0.01 
GeV^ < < 9.8 GeV^. In a particular analysis we shall 
also use the pp data at y/s = 27.5 GeV, in the region 5.5 
GeV^ < q^ < 14.2 GeV^. Some treatment of cosmic-ray 
information on pp total cross sections at y/s = 6-40 TeV 
is also presented. 

In Figure 1 it is displayed the experimental information 
available on pp and pp total cross sections from acceler- 
ators and cosmic-ray experiments. From that plot, it is 
clear that the mathematical description of the increase of 
the total cross sections at the highest energies is an open 
problem. As we shall discuss, the study of the effects of 
the discrepant points at the highest energies is one of our 
goals. Figure 2 shows the typical diffractive pattern that 
characterizes the differential cross section. We note that 
the data cover the region corresponding to 10 decades. In 
Figure 3 it is displayed the slope parameter from pp and 
pp scattering as function of the energy and determined in 
the region of small momentum transfer. In what follows 
we shall refer to these three figures as indicative of the 
empirical behavior of the quantities involved. 
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FIG. 1: Total cross sections on pp and pp from accelerator 
and cosmic-ray experiments (for complete list of references 
and tables see [lo|'). 



FIG. 3: The slope parameter as function of the energy and 
determined in the interval 0.01 < Itl < 0.20 GeV^. 



particle and particle-antiparticle interactions) . Analytic- 
ity and crossing allow the connections between real and 
imaginary parts of the scattering amplitude by means of 
dispersion relations. 

Several rigorous theorems and bounds may be deduced 
from the basic Principles and axiomatic quantum field 
theory. Among them, the Froissart-Martin bound con- 
cerns the increase of the total cross section stating that 
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oo. 



So 



(5) 



The Pomeranchuk Theorem treats the difference 
between cross sections for particle-particle (ab) and 
particle-antiparticle scattering (ab). The original form 
was deduced when it was believed that the cross section 
decreased to a constant value, and in this case cr^^j = ct^qj 
as s ^ oo. After the discovery of the rising of the cross 
section, Grunberg and Truong obtained the generalized 
or revised form of the Pomeranchuk Theorem, stating 
that 



'tot 
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FIG. 2: Differential cross section data and the diffractive pat- 
tern from pp elastic scattering at = 52.8 GeV. 



and this means that, 
reached, then 



if the Froissart-Martin bound is 



B. Principles, theorems and high-energy bounds 

For our purposes, we recall some principles and the- 
orems from axiomatic quantum field theories [ij. The 
basic Principles are: Lorentz Invariance, Unitarity (re- 
lated with the conservation of probability), Analyticity 
(related to causality) and Crossing (connecting particle- 



Aa ^ at - ^tt < C ''°°^f°* < Clogs. (6) 

By expressing the cross sections in terms of crossing 
even (-|-) and odd (— ) contributions. 
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we have I Act I 



\ ab _ ab | 



^ab 
'tot 



Hot - -toil = 2(T_ . Therefore, Act = 
(^tot - ^tot if and only if a_ 0. This possible 
odd contribution is named Odderon and the case of even 
dominance at asymptotic energies is associated with the 
Pomeron. 



C. Basic pictures 

Nearly all the phenomenological models, able to de- 
scribe the experimental data on elastic hadron scatter- 
ing, are based on the Optical/Geometrical Picture (s- 
channel) and/or the Exchange Picture (t-channel). The 
corresponding formulas may be obtained from the Partial 
Waves representation of the scattering amplitude, 



1=0 



where 5i is the phase shift. In what follows, we outline 
the main steps and formulas in both pictures. 



1. Optical/ Geometrical Picture 

From the partial wave representation, one considers the 
high-energy limit and the semi-classical approximation, 
so that the discrete angular momentum I may be replaced 
by the continuum impact parameter 6, 

l = kb-\. 

In turn, the discrete phase shifts 5i are replaced by the 
continuum eikonal function of h and s, x(s, b) and 



db... 



The scattering amplitude in this Eikonal Representa- 
tion, with azimuthal symmetry assumed, reads 



/■OO 

F{s, q)=ik I bdbJo{qb)[l - e 
Jo 

The quantity 

1 _ e'x(«.f') = r(s, b) 



_ pix(s,b)-\ 



(7) 



(8) 



is named Profile function. Prom Unitarity this function 
is related to the probability that an inelastic event takes 
place at b and s, the Inelastic Overlap function: 



Ginei{s,b) = \T{s,b)\^ -2Rer{s,b). 



(9) 



Since in the Eikonal representation 



Ginelis,b) = 1 



-2Imxis,b) 



(10) 



for Im x(s, 6) > we have Ginei{s, b) < 1, which implies 
in an automatically unitarized representation. 



2. Exchange Picture 

In this picture, from the partial wave representation, 
one considers the analytic continuation of the amplitude 
to complex angular momentum. In the asymptotic limit 
(s oo) and with symmetry connecting the crossed 
channels one arrives at the Watson-Sommerfeld-Gribov- 
Reggc representation for the scattering amplitude, ex- 
pressed as a sum over the poles of the amplitude (the 
Regge poles), as outlined in what follows. 

As it is known at high energies the number of par- 
tial waves is large, and one way to circumvent that is to 
transform the sum of partial waves into a complex inte- 
gral, and then use the residues theorem to obtain a new 
sum, but involving only the number of residues: 



oo „ N 



1=0 



Detailed calculation allows one to obtain the following 
representation for the scattering amplitude, 



N 



1=1 



sin iTai{k) 



where BI{k, 6) is called the Background integral. By 
considering the high-energy limit (then BI 0) and 
crossing (exchange four-momenta p ^ 4^ ^ —p) we 
can replace the crossing channel variable {9 s) 



cos 6=1 — 



2s 



also. 



Aw? - 1 



2'r{i + 1/2) 



oc —s as s — > oo, 



V^r{i + 1) 



X for X — > ±oo, 



and grouping all the s-independent quantities in a func- 
tion )C{t) we have 



-Pa(t) (- cos 61) = /C(t)s"(*) for s ^ oo. 

Rearranging the terms we arrive at a descending 
asymptotic series in powers of s, with leading contribu- 
tion: 
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FIG. 4: The Chew-Frautschi plot for some mesons and reso- 
nances. 



F(s,t)=7(0^(t)s"(*), (11) 

where 7(t) is the residue function, ^(t) the signature fac- 
tor and a(t) = a(0) -I- act the trajectory function. This 
last function connects the spin and masses through the 
Chew-Frautschi plot, as exemplified in Fig. 4. In this pic- 
ture the interaction of the colliding particles is basically 
interpreted in terms of exchanges of Regge poles (also 
Regge cuts) and the Pomeron (an ad hoc trajectory with 
intercept nearly above 1). We note that, as constructed, 
the exchange picture is intended for asymptotic energies. 



At high energies, the standard singly subtracted inte- 
gral dispersion relations, with poles removed, are given 
by 

2„2 /■ + 00 . 

RcF+{s)=K + —P ds'-—^ ^ImF+(s') 

J so S (S - S ) 

(13) 

and 

2s 1 
ReF.{s) = -P ds'— -h^^F^is'), (14) 

I" J so (S - S ) 

where K is the subtraction constant and, for pp and pp 
scattering, sq — '--^ 1.8 GeV^. 

In this section we review some results obtained through 
this approach. We start with the replacement of the 
above integral forms by derivative operators (Deriva- 
tive Dispersion Relations) and then we discuss the use 
of analytic models (Reggeons, Pomeron, Odderon) for 
parametrizations involving the total cross section and the 
p parameter, the determination of bounds for the soft 
Pomeron intercept, and the practical role of the subtrac- 
tion constant. 

In what follows we are mainly concerned with the pp 
and pp elastic scattering, since for particle and antipar- 
ticle interactions they correspond to the highest energy 
interval with available data and are the only set includ- 
ing the cosmic-ray information on total cross sections {pp 
scattering) . As commented before, the experimental data 
available on the total cross sections (Figure 1) are char- 
acterized by discrepant experimental information at the 
highest energies, and one of our aims is to investigate the 
effects of these discrepancies in the context of the ana- 
lytic models. This concern permeates all the discussion 
in this Section. 



III. ANALYTIC APPROACH 

The Analytic Approach for elastic hadron-hadron scat- 
tering is based on general principles and theorems from 
Quantum Field Theory. It is characterized by analytical 
parametrizations for the imaginary part of the forward 
amplitude, together with the use of dispersion relation 
techniques. The central point is the simultaneous inves- 
tigation of the total cross section (imaginary part of the 
scattering amplitude, Eq. (2)) and the p parameter (con- 
nected with the real part of the amplitude, Eq. (3)). 

For particle-particle and particle-antiparticle interac- 
tions, dispersion relations are consequences of the prin- 
ciples of Analyticity and Crossing. In this context, they 
correlate real and imaginary parts of crossing even (-I-) 
and odd (— ) amplitudes, which in turn are expressed in 
terms of the scattering amplitudes for a given process 
and its crossed channel, for example, a -I- 6 and a + h: 

Fab = F++F^, F,-, = F+-F^. (12) 



A. Derivative Dispertion Relations 

The use of dispersion relations in the investigation of 
scattering amplitudes may be traced back to the end of 
fifties, when they were introduced in the form of Inte- 
gral Dispersion Relations (IDR). Despite the important 
results that have been obtained since then, one limita- 
tion of the integral forms is their non-local character: in 
order to obtain the real part of the amplitude, the imag- 
inary part must be known for all values of the energy. 
Moreover, the class of functions that allows analytical 
integration is limited. 

In the last years, we have investigated the applicabil- 
ity of Derivative Dispersion Relations (DDR) in place 
of integral forms |ip, 12, 13_, J^, 16J. In Reference 
|16| we present a recent review on different results and 
statements related to this replacement, and a discussion 
connecting these different aspects with the corresponding 
assumptions and classes of functions considered in each 
case. 
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In particular, we have shown that for the class of func- 
tions which are entire in the logarithm of the energy (as 
is the case of analytic models at high energies) it is pos- 
sible to expand the integrand in the above formulas and 
by considering a high-energy approximation, represented 
by So = 2m^ 0, to integrate term by term. In that 
case, as demonstrated in detail in the derivative dis- 
persion relations with one subtraction reads 



ReF+(s) K 

= h tan 

s s 



2 dins 



ReF_(s) 



= tan 



2 V ' dins 



ImF+(s) 



ImF_(s) 



(15) 



(16) 



where the series expansion is implicit in the tangent op- 
erator. From this deduction one arrives to three formal 
results: (1) the subtraction constant is preserved when 
the IDR are replaced by DDR and, therefore, in princi- 
ple, can not be disregarded in fit procedures; (2) except 
for the subtraction constant, the DDR with entire func- 
tions in the logarithm of the energy do not depend on any 
additional free parameter; (3) the only approximation in- 
volved in the replacement concerns the lower limit in the 
IDR (13-14), namely sq = 2m^ 0, which represents a 
high-energy approximation. In the next two subsections 
we discuss some uses of the DDR with analytical models, 
and in the third subsection we return to the replacement 
of IDR by DDR, investigating the important role of the 
subtraction constant from a practical point of view. 



Taken separately these two sets suggest different scenar- 
ios for the increase of the total cross section, as previously 
discussed in [ll[ll[l3. 

Based on these considerations, it is important to in- 
vestigate the behavior of the total cross section by tak- 
ing into account the discrepancies that characterize the 
cosmic ray information. To this end, in 'lOl| we have 
considered two ensembles of data and experimental in- 
formation, as follows: 

• Ensemble I: pp and pp accelerator data + Akeno -t- 
Fly's Eye; 

• Ensemble II: pp and pp accelerator data -I- Nikolaev 
+ GSY. 

To some extent, ensemble I represents a kind of high- 
energy standard picture and ensemble II a nonstandard 



2. Analytic Models 

With analytical parametrizations for pp/pp total cross 
sections, the connections with the p parameter, Eq. (3), 
arc obtained by defining the associated crossing even and 
odd quantities, 



cr±(s) 



-'tot "tot 



(17) 



using the high-energy normalization for the Optical The- 
orem, 



B. Basic Models 

In this Subsection we make use of two basic and well 
known parametrizations for the total cross sections and 
investigate the effects of the discrepancies in the experi- 
mental information from cosmic-ray experiments. 



1. Ensembles 

In the cosmic-ray region, 6 TeV < -y/s < 40 TeV, the 
discrepancies on the total cross section information are 
due to both experimental and theoretical uncertainties 
in the determination of crf^j from p-air cross sections. 
The situation has been recently reviewed in detail in [lOf , 
where a complete list of references, numerical tables and 
discussions are presented. 

From Fig. 1 we see that, despite the large error bars in 
the cosmic-ray region, we can identify two distinct sets 
of estimations: one corresponding to the results by the 
Fly's Eye Collaboration (Fly's Eye) together with those 
by the Akeno Collaboration (Akeno); the other set asso- 
ciated with the results by Gaisser, Sukhatme, and Yodh 
(GSY) together with with those by Nikolaev (Nikolaev). 



(Ttot{s) 



ImF(s,0) 



(18) 



and the DDR given by Eqs. (15) and (16). 

In |ipj we have considered two different parametriza- 
tions for the total cross sections, one introduced by Don- 
nachie and Landshoff 19] and other by Kang and Nico- 
lescu [201. The main difference concerns the asymptotic 
limits, which allow the dominance of an even amplitude 
(Pomeron) or the odd amplitude (Odderon), respectively. 
In this way, we may contrast these possibilities with the 
standard and non-standard pictures represented by En- 
sembles I and II. 

The Donnachie- Landshoff (DL) parametrization for 
the total cross sections is expressed by 



aZ{s)=Xs^ + Zs-^, (19) 



where the first contribution is associated with a single 
Pomeron exchange (universal) and the second one with 
Reggeon exchange. With the procedure explained above, 
we obtain the analytical connections with the p parame- 
ter for pp and pp scattering: 
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K 



K 



X tan 



V 2 



p''^(«(.) = - + 



X tan 



2 V 2 y 2 V 2 



small, so that we may replace crfoj « crfoj 
asymptotically, 



CTtot(s), then. 



T^CTtotis) 



{A^lns + ABln^s} . 



This means that, depending on the fit results, there may 
be a change of sign in Ap, with p^^ becoming greater than 
at some finite energy. Therefore, the case of a crossing 
either in atot or p is a sign of the odderon contribution in 
the imaginary or real part of the amplitude, respectively. 



From the above formulas, since ?7 > 0, this model pre- 
dicts that, asymptotically (s oo), 



3. Fits and Results 



The parametrization for the total cross sections intro- 
duced by Rang and Nicolescu (KN), under the hypothesis 
of the Odderon, is given by 



aZ{s) = Ai + Bi\ns + k In s, 



aZis) = A2 + B2lns + klii\s 
and the connections with p read 



2R 



pPP{s)aZ^{s) = - + - 
s Z 



TTK H in s - 



2 

B2 — Bi 
2tt 



1 2 2i? 

In s 



1/2 ' 



ffV[s)aZt{s) = - + - 
nk 



K TT / Bi+B: 



s 2 
In s — 



A2-Ai \ , f B2 - Bi N 2 
2tt ^ 



Differently from the previous case, this model predicts 
that the difference between the two cross sections is given 

by 



Act ={A2-Ai) + {B2 -Bi)lns + 2Rs-^''^ 
AA + A_Blns (asymptotically), 

so that, if AA ^ and/or AB ^ 0, the total cross section 
difference may increase and aZt may even become greater 
than crZt, depending on the values and signs of AA and 
Ai?, which is formally in agreement with the theorems 
of Sec. II. B. Moreover, if AA and AB are sufficiently 



We have performed 16 different fits through the pro- 
gram CERN-MINUIT. In these fits we have used both 
ensembles I and II and both the DL and KN models. For 
each of these four possibilities we have performed global 
and individual fits to atot and p and, in each case, we 
either considered the subtraction constant _?C as a free fit 
parameter, or assumed K = 0. 

All the results are presented and discussed in detail in 
Ref. 0^'' main conclusions are the following: (1) 

Despite the small influence from different cosmic-ray es- 
timations, the results allow to extract an upper bound 
for the soft Pomeron intercept: 1 + e = 1.094; (2) al- 
though global fits present good statistical results, in gen- 
eral, this procedure constraints the rise of atot', (3) the 
subtraction constant as a free parameter affects the fit 
results at both low and high energies; (4) independently 
of the cosmic-ray information used and the subtraction 
constant, global fits with the Odderon parametrization 
predict that, above ^/s « 70 GeV, Ppp{s) becomes greater 
than Ppp(s), and this result is in complete agreement with 
all the data presently available. That result is displayed 
in Fig. 5 and we can infer ppp = 0.134 ± 0.005 at 
^ = 200 GeV and 0.151 ± 0.007 at 500 GeV (BNL 
RHIC energies). 



C. Non-degenerate Meson Trajectories 

The DL parametrization referred above assumes de- 
generacies between the secondary reggeons, imposing a 
common intercept for the C = -1-1 (02,72) and the 
C = ~1 {ijJ,p) trajectories (see Fig. 4). More recently, 
analysis treating global fits to atot and p have indicated 
that the best results are obtained with non-degenerate 
meson trajectories. In this case the forward scattering 
amplitude is decomposed into three reggeon exchanges, 
F{s) = F3>(s) + F^^/f^{s) + TF^,p{s), where the first 
term represents the exchange of a single soft Pomeron, 
the other two the secondary Reggeons and r = -1-1 (—1) 
for {pp) amplitudes. Using the notation Q;p(0) = 1 + e, 
Q!_|_(0) = 1 — 77+ and Q;_(0) = 1 — 7]_ for the intercepts of 
the Pomeron and the C — +1 and C — ~\ trajectories. 
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FIG. 5: Simultaneous fits to atot{s) and p{s) through the KN 
parametrization with K — and ensembles I (dotted curves 
for pp and dashed for pp) and II (solid curves for pp and dot- 
dashed for pp) p^. 



respectively, the total cross sections, Eq. (18), for pp and 
pp interactions are written as 



atot{s) = Xs' + Y+ 5-"+ + tF_ s-"- (20) 

and the connection with the p parameter by means of 
DDR is similar to that displayed in the last subsection. 

Making use of this parametrization, in this section 
we present the determination of extrema bounds for the 
Pomeron intercept Il5!| and a practical analysis on the 
replacement of IDR by DDR together with a discussion 
on the role of the subtraction constant ligl. 



1. Extrema Bounds for the Pomeron Intercept 

In order to analyze the extrema effects in the soft 
Pomeron intercept due to discrepancies in the experimen- 
tal data, we performed a detailed analysis including the 
highest and the lowest values of the total cross section 
from both accelerators and cosmic-ray experiments. 

As it is well known, in the accelerator region, the con- 
flict concerns the results for crf^j at =1.8 TeV re- 
ported by the CDF Collaboration and those reported by 
the E710 and the E811 Collaborations (Fig. 1). In the 
cosmic-ray region, as we have discussed, the highest pre- 
dictions for cr^Qj concern the result by Gaisser, Sukhatme, 
and Yodh together with those by Nikolaev. In order 
to treat the lowest estimations in the cosmic-ray region, 
we consider the results obtained by Block, Halzen, and 
Stanev (BHS), by means of a QCD-inspired model. As 
discussed in |Tq|, the reason for this choice is that, al- 
though the extracted crf^^^s) shows agreement with the 
Akeno results, it is about 17 mb below the Fly's Eye value 
at 30 TeV and therefore may be considered as a extreme 
lower estimate. All the numerical tables and references 
can be found in 0| . 

In this case we have considered the following ensem- 
bles of experimental information. First we only consider 
accelerator data in two ensembles with the following no- 
tation: 

• Ensemble I: af^^ and af^t data (10 < < 900 GeV) 
-I- CDF datum (^/i = 1.8 TeV); 

• Ensemble II : af^^ and af^t data (10 < ^/i < 900 
GeV) + E710/E811 data (^i = 1.8 TeV). 

Ensemble I represents the faster increase scenario for 
the rise of atot from accelerator data and ensemble II 
the slowest one. These ensembles are then combined 
with the highest and lowest estimations for af^^ from 
cosmic-ray experiments, namely, the Nikolaev and the 
Gaisser, Sukhatme, and Yodh (NGSY) results and the 
Block, Halzen, and Stanev (BHS) results, respectively. 
These new ensembles are denoted by 

• I -I- NGSY 

• II + BHS 

As in the previous analysis, we have considered both 
individual fits to atot, and simultaneous fits to atot and 
p, either in the case where the subtraction constant is 
considered as a free fit parameter or assuming K = 0. 

From this analysis, in the case of only accelerator data, 
we could infer the following upper and lower values for 
the Pomeron intercept: aff(O) = 1.098±0. 004 (global fits 
to ensemble I, with K = 0) and ^^(O) = 1.085±0.004 (in- 
dividual fit to atot from ensemble II), with bounds 1.102 
and 1.081, respectively. Adding the cosmic-ray informa- 
tion, we inferred the following upper and lower values: 
aff(O) = 1.104 ± 0.005 (individual fit to atot from en- 



9 




FIG. 6: Fastest and slowest increase scenarios for the rise 
of the total cross section through parametrization (20) and 
allowed by the experimental information available: fits to en- 
sembles I + NGSY (solid) and II (dashed) [IJ. 



semble I + NGSY) and ai.(0) = 1.085 ± 0.003 (global 
fits to ensemble II + BHS and if as a free fit parameter 
or individual fit to atot from this ensemble) , with bounds 
1.109 and 1.082, respectively. Therefore we may infer the 
following extrema bounds for the soft Pomeron intercept: 



^.pper^Q^ = 1.109, 



^(0) = 1.081. 



2. IDR, DDR and the Subtraction Constant 

As commented before, we have shown in Ref. [T^ 
that for entire functions in the logarithm of the energy 
the only approximation involved in the replacement of 
IDR by DDR concerns the lower limit sq in the IDR: 
the high-energy condition is reached by assuming that 
So = Irn? ^ in Eqs. (13-14). In that paper we have 
investigated the practical applicability of the DDR and 
IDR in the context of the Pomeron-reggeon parametriza- 
tions, with both degenerate and non-degenerate higher 
meson trajectories. By means of global fits to atot{s) 
and p{s) data from pp and pp scattering, we have tested 
all the 16 important variants that could affect the fit re- 
sults, namely the number of secondary reggeons, energy 
cutoff (5 and 10 GeV), effects of the high-energy approx- 
imation connected with the subtraction constant and the 
analytic approach using both DDR and IDR with fixed 
So- Our results led to the conclusion that the high-energy 
approximation and the subtraction constant affect the fit 
results at both low and high energies. This effect is a con- 
sequence of the fit procedure, associated with the strong 
correlation among the free parameters. 

A striking novel result concerns the practical role of 
the subtraction constant. We have shown that, with the 
Pomeron-reggeon parametrizations, once the subtraction 
constant is used as a free fit parameter, the results ob- 
tained with the DDR and with the IDR (with finite lower 
limit. So — 2m^) are the same up to 3 significant figures in 
the fit parameters and /DOF. This conclusion, as we 
have shown, is independent of the number of secondary 
reggeons (DL or extended parametrization) or the energy 
cutoff {^/s = 5 or 10 GeV). In Table I we display the fit 
results with the extended parametrization and cutoff at 
10 GeV. 



Figure 6 shows the total cross sections with parametriza- 
tion (20) and the above extrema bounds, together with 
the experimental information available. 

Extensions of these extrema hounds for the pomeron 
intercept to meson-p, gamma-p and gamma-gamma scat- 
tering have been discussed in (2l|. By means of global 
fits to total cross section data it is shown that these 
bounds are in agreement with the bulk of experimen- 
tal data presently available, and extrapolations to higher 
energies indicate different behaviors for the rise of the 
total cross sections. 

We have also obtained new constrained hounds for 
the Pomeron intercept from spectroscopy data (Chew- 
Frautschi plots) and have extended the analysis to 
baryon-p, meson-p, baryon-n, meson-n, gamma-p and 
gamma-gamma scattering |22|. It is also presented tests 
on factorization and quark counting rules with both ex- 
trema and constrained bounds (asymptotic energy re- 
gion). In particular, at 14 TeV (CERN LHC) the extrema 
and constrained bounds allow to infer atot = 114 ±25 mb 
and 105 ± 10 mb, respectively. [2^ . 



TABLE I: Simultaneous fits to atot and p through the ex- 
tended parametrization, V^min ~ 10 GeV (154 data points), 
with K as. a. free parameter and using IDR with lower limit 
So = 2m^ and DDR [l^. 





IDR with So = 2m^ 


DDR 


X 


19.57 ± 0.79 


19.58 ± 0.78 


y+ 


66.0 ± 6.7 


66.0 ± 6.6 


Y- 


-29.2 ± 4.0 


-29.2 ± 4.0 


t 


0.0897 ± 0.0033 


0.0897 ± 0.0033 


V+ 


0.380 ± 0.033 


0.380 ± 0.033 


V- 


0.520 ± 0.025 


0.520 ± 0.024 


K 


-14 ± 48 


104 ± 58 


xVdof 


1.10 


1.10 



IV. MODEL INDEPENDENT ANALYSIS 

This kind of analysis is characterized by model inde- 
pendent parametrizations of the experimental data in- 
volved and the extraction of empirical or semi-empirical 
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information that can contribute with the development of 
phenomenological models and the underlying theory. In 
this section we review some results we have obtained in 
the investigation of pp and pp differential cross section 
data (unconstrained and constrained fits, as will be ex- 
plained) and the correlations between the experimental 
data on total cross section and the slope parameter (Eq. 
(4)). 

A. Differential Cross Section 

Several authors have investigated elastic hadron scat- 
tering by means of parametrizations for the scattering 
amplitude and fits to the differential cross section data, 
Eq. (1). The extraction of the Profile, Eikonal and 
Inelastic Overlap functions in the 6-space (impact pa- 
rameter) and, in some special cases, the Eikonal in the 
q^-space, has led to important and novel results related 
with geometrical aspects (radius, central opacity), differ- 
ences between charge distributions and hadronic matter 
distributions, existence or not of eikonal zeros in the q^- 
space and, more recently, connections with pomerons, 
reggeons and nonperturbative QCD aspects. In Ref. 
we present a review and a critical discussion on the main 
results concerning this kind of analysis and also a wide 
list of references to outstanding works. 

The basic input in all these analyses is the parametriza- 
tion of the scattering amplitude as a sum of exponentials 
in (as empirically suggested by the diffractive pattern 
shown in Fig. 2) and fits to the differential cross section 
data. This parametrization allows analytical expressions 
for the Fourier transform of the amplitude, providing also 
analytical expressions for the quantities of interest in the 
6-space. 

In the next two subsections we review the results we 
have obtained by means of unconstrained fits (fit param- 
eters completely free, without extracted dependences on 
the energy) |23, UM i ^"^d discuss some research in course 
related to constrained fits (including dependences on the 
energy which are based on empirical information) |l27l| . 

1. Unconstrained Fits and the Eikonal 

In the high energy region, ^/s >10 GeV, differential 
cross section data are available at ^/s — 13.8, 19.5, 23.5, 
30.7, 44.7, 52.8 and 62.5 GeV for pp scattering and at 
y/s = 13.8, 19.4, 31, 53, 62, 546 and 1800 GeV for pp 
scattering. Data from pp scattering also exists at ^/s — 
27.5 GeV and 5.5 < < 14.2 GeV^ (but not on atot and 
p), and as we shall show, that set plays a fundamental role 
in our analyse. See 23| for a complete list of references. 

As discussed in 23] two main problems are typical of 
model independent analysis of the differential cross sec- 
tions: 

(1) Experimental data are available only over finite 
regions of the momentum transfer (which in general are 



small, q"^ < 7 GeV^) and the Fourier transform demands 
integration from = to infinity. This means that 
any fit is biased by extrapolations and although some 
extrapolated curves may look unphysical, they can not 
be excluded on mathematical grounds. 

(2) The exponential parametrization allows analytical 
determination of the quantities in the 6-space (profile, 
inelastic, eikonal functions) and also the statistical un- 
certainties, by means of error propagation from the fit 
parameters. However, in this case, the translation of 
the eikonal from b-space to the g^-space can not be an- 
alytically performed and neither the error propagation 
(through standard procedures). As a consequence, the 
unavoidable uncertainties from the fit extrapolations can 
not, in principle, be taken into account. 

In what follows we review a model independent ap- 
proach able to minimize the above two problems. 

- Fit Procedure 

In order to trea.t problem (1) we have used the fol- 
lowing procedure Since it is known that for large 
t the experimental data do not depend on the energy at 
13.8 GeV < y/s < 62 GeV and that there exist data at 
= 27.5 GeV in the region 5.5 GeV^ < < 14.2 
GeV^, we have selected two ensembles of pp and pp dif- 
ferential cross section data: 

• Ensemble I: experimental data at each energy; 

• Ensemble II: Ensemble I -I- data at ^/s = 27.5 GeV. 
For the scattering amplitude we have introduced the 

following model independent analytical parametrization 
for both real and imaginary parts: 

2 n 

F{s,q)^{fiY.a,e'f'^'^'} + t{Y,a,e-^^'''}, (21) 

,= ^^±a,. (22) 
ai + a2 ^ 

With the experimental p value at each energy the fits 
to the differential cross section data have been performed 
through the CERN-MINUIT routine and the vahdity or 
not of ensemble II is checked by means of the MINUIT 
output and standard statistical interpretation of the fit 
results (DOF, confidence levels). 

For pp scattering we have found that the data at ^/s — 
13.8 GeV are not compatible with ensemble II. In the 
case of pp scattering none of the data sets are compatible 
with ensemble II. Therefore, in what follows, ensemble II 
(data at y/s — 27.5 GeV added) corresponds only to pp 
scattering at 6 energies: 19.5, 23.5, 30.7, 44.7, 52.8 and 
62.5 GeV. 

From the error matrix (variances and covariances), 
/DOF and confidence intervals, we infer the best val- 
ues for the parameters and corresponding errors Aa^, 
A/3j. By means of standard error propagation, the un- 
certainties in the free parameters, Aa^, A/3j, (j = 1, 2, 
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FIG. 7: Regions of uncertainties (limited by the solid lines) in 
fits to pp differential cross section data, Eq. (23 ], a t -^s = 23.5 
GeV with ensembles I (below) and II (above) |2^ . 



FIG. 8: Typical parametrization for the generated remainder 
R{s, b) by means of Eq. (26) 



...) have been propagated to the scattering amplitude, 
and then to the differential cross section, providing 



da 



± A 



da 
d^ 



(23) 



By adding and subtracting the corresponding uncertain- 
ties we may estimate the confidence region associated 
with all the extrapolations, which cannot be excluded on 
statistical grounds. A typical result with ensembles I and 
II is illustrated in Fig. 7, for pp scattering at \/s= 23.5 
GeV. We see that, as expected, the effect of adding the 
experimental data at ^/s= 27.5 GeV (when statistically 
justified) is to reduce drastically the uncertainty region. 
That result will be fundamental in the extraction of the 
empirical information on the eikonal, as shown in what 
follows. 

- Eikonal in the momentum transfer space 
By means of the Fourier transform, Eqs.(7-8), the 
parametrization (21-22) provides analytical expressions 
for the real and imaginary parts of the Profile function, 
Tfi{s^h) and r/(s,5), and also the associated uncertain- 
ties. From the fit results, together with error propaga- 
tion, we have found that 



«1, 



and therefore, the imaginary part of the eikonal may be 
approximated by 



X/(s,5) w In 



1 



l-TR{s,h) 



(24) 



and the uncertainty Ax/ determined directly from A^i^ 
through propagation. 

The next step is to go to the momentum transfer space 
and concerns problem (2): the Fourier transform can not 
be performed analytically and therefore also the error 
propagation. For this reason we used a semi-analytical 
method as follows. Expanding the above equation, we 
express the remainder of the series as 



i?(s, 6) = ln[ 



1 



TR{s,b) 



(25) 



and then fit the numerical points (MINUIT) by a sum of 
Gaussians in the impact parameter space: 



i=i 



(26) 



A typical result is displayed in Fig. 8. 

With this, the errors AAj and ABj may be propagated 
determining Ai?/it(s, g) and then X7(s, g)±Ax/(s, g). In 
order to check the results and approximations, we per- 
formed also numerical integration through the NAG rou- 
tine. 

- Results 

One of the main results extracted from this analysis is 
the statistical evidence of eikonal zeros in the momentum 
transfer space, first presented in . In order to investi- 
gate the position of the zeros and, mainly, to determine 
the uncertainties in its values, we consider the expected 



behavior of xi large q , namely x/ 



In Fig. 



(9) we show a typical plot of the quantity q^Xi{s,q) as 
function of . The shaded areas correspond to the un- 
certainties obtained from error propagation. This exam- 
ple shows clearly the role and the effect of data at large 
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FIG. 9: The eikonal in the transfer momentum space (multi- 
phed by t*) for pp at ^/s = 30.7 GeV with ensemble I (above) 
and II (below) |2^ . 



values of the momentum transfer. In fact, within the 
uncertainties, ensemble I does not allow to infer a zero, 
but with ensemble II, we find statistical evidence for the 
change of sign. 

From plots like that we can determine the positions 
of the zeros and the associated errors from the extrems 
of the uncertainty region (in general not symmetrical). 
The position of the zero can also be obtained from the 
numerical method, but without uncertainties. 

In Figure 10 it is shown the position of the zeros as 
function of the energy determined by means of both 
the semi-analytical (with uncertainties) and numerical 
(without uncertainties) methods. Despite the system- 
atic difference on the values with these methods, we may 
conclude that the position of the zero decreases as the 
energy increases. Roughly, : 8.5 —> 6.0 GeV^ as 

: 20 ^ 60 GeV. 
- Discussion 

As reviewed in '2^ , there has been previous indication 
of eikonal zeros in the momentum transfer space, but 
without associated uncertainties. Our first statistical ev- 
idence, published in 1997, indicated the position of the 
zero at = 7 ± 2 GeV^ In 2000, experiments per- 
formed at the Jefferson Laboratory, on electron-proton 
scattering, have indicated an unexpected decrease of the 



> 
o 



- - o semi-analytical 
X X numerical 



20 30 40 50 60 70 

s'" (GeV) 

FIG. 10: Position of the eikonal zero in the momentum trans- 
fer space as function of the energy |2^. 



ratio between the electric and magnetic proton form fac- 
tors as the momentum transfer increases from 0.5 to 5.6 
GeV^. Moreover, extrapolations from empirical fits in- 
dicate a change of sign (zero) in this ratio, just at « 
7.7 GeV^. Since for pp/pp scattering the eikonal is con- 
nected with the hadronic matter form factor (see Sec. 
V, Eq. (38)), the above results on the position of the 
zeros suggest novel and important insights on possible 
correlations between hadronic and electromagnetic inter- 
actions. We discuss that subject in calling attention 
to the possibility of hadronic form factors depending on 
the energy. 

We have also obtained the value of the imaginary part 
of the Eikonal at zero momentum transfer, that is, the 
central opacity. The results are displayed in Fig. 11. 
As discussed in one naive way to test these results 
is with the Glauber model for the scattering involving 
hadrons A and B, and the Optical Theorem at the ele- 
mentary level. In that case, the elementary cross section 
may be expressed by 



47r 



NaNb 



Xi{s,q^ 0), 



where Na and Nb are the number of constituents in 
hadrons A and B, respectively. If we take Na = Nb = 3 
we obtain aeiem ~ 6 mb at the ISR region, a result in 
agreement with other estimations. 

In Ref. I23 we discuss the applicability of our re- 
sults in the phenomenological context, outlining some 
connections with nonperturbative QCD and presenting 
a critical review on the main results concerning "model- 
independent" analyses. 
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FIG. 11; Imaginary part of the eikonal at zero momentum 
transfer, as function of the energy, from analysis on pp and 
PP scattering (23(1 . 



Y,c^^{s)eM-P^{s)q% (29) 



4=1 



and, based on the above arguments, we introduce the 
following general dependences on the energy 



ai{s) ^ Ai + Biln{s) + C,ln{s) 



(30) 



for pp scattering and 



a^{s) = A, + B,\n{s) + ain\s) 
A(s) = A: + -B,ln(s) 



(31) 



for pp scattering, where i — 1,2, ...n. Defining the cross- 
ing even (+) and odd (-) amplitudes. 



2. Constrained Fits and Energy Dependence 

Despite the results obtained with the parametrization 
discussed in the last subsection, due to the fit procedure, 
we do not have the dependence on the energy of the free 
parameters at, Pi. Presently, we are investigating that 
subject and we review here some preliminary results |27j . 

The energy dependence has been introduced according 
to some empirical information: the increase of the total 
cross section and of the slope parameter with In^ s and 
Ins, respectively (see Figs 1 and 3). Let us consider the 
standard exponential parametrization for the imaginary 
part of the amplitude, normalized as 



lmF{s,q^) 



i=l 



ai exp[-/3ig^ 



(27) 



At = 0, from the optical theorem, Eq. (18), we 
expect a dependence of the parameters with In^ s, 
and the slope represented by the parameters (3i with In s. 
These are the central choices in our approach. In order 
to treat pp and pp scatterings, in agreement with An- 
alyticity and Crossing, we introduce crossing even and 
odd amplitudes and make use of the derivative disper- 
sion relations, Eqs. (15) and (16), to connect real and 
imaginary parts of the amplitudes involved. 

Specifically, for the imaginary part of the scattering 
amplitude we consider the parametrizations 



Im Fpp(s,g2) 



= y]"j(s) exp[-A(s)(j^ 



=1 



(28) 



lmF^{s,q'^) 



Im Fpp(s,g2) - Im i^pp(s,g^) 



(33) 



the corresponding real parts can be determined by means 
of the leading terms of the DDR, Eqs. (15-16), and so 
the corresponding real parts of the pp and pp amplitudes. 
With these analytic amplitudes we obtain the differential 
cross section: 



da 



167rs^ 



\ReF{s,q^) +ilTJiF{s,q 



2\|2 



(34) 



In order to treat simultaneous fits to pp and pp data 
we have considered only sets available at nearly the same 
energy, namely ^/s ~ 19.5, 31, 53 and 62 GeV. As a 
preliminary test we make use of data at the diffraction 
peak, outside the Coulomb-nuclear interference region, 
0.01 GeV^ <q^ < 0.5 GeV^, and the data providing the 
optical point. 



da(s,g2 = 0) cTtot{l + P^) 



dq^ 



(35) 



We have performed simultaneous fits to the experimen- 
tal data through the MINUIT program. For this ensem- 
ble we used only two exponentials for the imaginary part 
of the amplitude, obtaining good reproduction of all the 
data analyzed, as shown in Fig. 12. In Ref. |^ we also 
display the predictions for the differential cross sections 
at the RHIC energies. We are presently investigating the 
extension of the analysis to the region of higher momen- 
tum transfer. 
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FIG. 12; Differential cross section at the diffraction peak: fit 
results and experimental data (displaced by a factor of 10) 



B. Total Cross Sections and Slopes 

Another important quantity that characterizes the 
elastic hadron-hadron scattering is the slope parameter, 
defined in Eq. (4). In practice it may be determined 
by means of fits to the hadronic differential cross section 
data in the region of small momentum transfer, with the 
parametrization 



da 



da 
'dt 



-B\t\ 



(36) 



and, in general, it is connected with p and atot through 
fits in the region of Coulomb interference 0. The slope 
and the total cross section are also important quanti- 
ties in the determination of crf^f from (cosmic-ray 



experiments) but, as commented before, the procedure is 
strongly model dependent. One reason is associated with 
the use of the Glauber multiple diffraction formalism, in 
which atot{s) and B[s) take part in the parametrization 
of the elastic amplitude. 



Fff(s,t)cxafP(s)exp 



B{s)t 



(37) 



As commented in [I^ and |23 , different models predict 
different relations between atot{s) and B{s) and that is 
mirrored in the final value of the cross section, contribut- 
ing to the discrepancies already discussed. 

Based on the above observations, we have investigated 
the possibility to extract an empirical correlation be- 
tween the experimental data on atot{s) and B{s), from 
pp and pp scatterings. 

For the slope parameter, we have selected the data 
above the region of Coulomb-nuclear interference and be- 
low the "break" in the hadronic slope at the diffraction 
peak (localized at |t| ~ 0.2 GeV^ at the ISR and CoUider 
energies), namely 0.01 < \t\ < 0.20 GeV^ (Fig. (3)). In 
this region, the differential cross section data are well 
fitted by a single exponential and therefore there is no 
change in the slope associated with the i-dependence. For 
each energy we have compiled the corresponding data on 
the total cross section. 

Once more, the choice for a parametrization was based 
on the empirical observation that at high energies B{s) 
increases with the logarithm of s. Since the Kang- 
Nicolescu parametrization for the total cross section is 
expressed in terms of Ins (Sec. III.B.2), we replaced this 
dependence by the slope parameter: 



ci + C2B + c^B^ 



(s) = ci + C2B + csB'^ 



046 



-B/2 



(38) 



where Ci, i — 1,2,3,4 are free fit parameters. That is a 
strictly mathematical choice, having nothing to do with 
the physics or model concept behind the Kang-Nicolescu 
parametrization. 

Fits to the experimental data have been performed 
with the CERN-MINUIT program and the results are 
displayed in Fig. 13. It is expected that extrapolations 
to cosmic-ray energies may be useful in the determina- 
tion of the pp total cross section from p-air cross section, 
allowing to connect atot — S in an almost model indepen- 
dent way. We are presently investigating this subject. 

In Ref. [25 we also made use of the Donnachie- 
Landshoff parametrization, which predicts a faster in- 
crease of the total cross section as function of the slope 
parameter. Moreover, in we also present a critical 
discussion on the recent measurement of the slope pa- 
rameter at the BNL RHIC, at 200 GeV, by the pp2pp 
Collaboration. We call attention to the fact that the 
combination B = 16.3 ± 1.8 GeV~^ and atot = 51.6 
mb, indicated by the pp2pp analysis, is in disagreement 
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FIG. 13; Total cross sections in terms of the slope, and the 
parametrization (38) '28'] ■ 



with the general trend for the behaviors of atot and B. 
If this "peer" is correct, new physics is necessary. Us- 
ing the above B value as input in our parametrizations, 
the corresponding values of the total cross sections show 
agreement with the atot versus B data. However, these 
inferred values for atot indicate new physics when plot- 
ted as function of the energy. We conclude that if this 
measurement is correct and represents an hadronic quan- 
tity, its high value may indicate a "break" in the slope 
near 0.02 GeV^, a phenomenon that was never observed 
in both pp and pp scattering, at ^/s < 62.5 GeV and 
^/s < 1.8 TeV, respectively and therefore, once more, 
new physics is necessary. 



- Elastic and Inelastic Channels 

Through the Unitarity and the Inelastic Overlap Func- 
tion, defined in Sec. II. C, we can connect elastic and 
inelastic scattering. This is done by expressing the topo- 
logical cross section for producting an even number n of 
charged particles at s in terms of Gin'. 



an{s) = J (fh an{s,b) ^ J 



cr„(s,5) 



Gin{s,b) 



If n(s) and < n >(,j) are the hadronic and averaged 
multiplicities, respectively, by introducting the KNO 
variable Z = n{s)/ < n >(s), the hadronic multiplicity 
distribution may be expressed by 



<l>(s, Z) =< n > 



anjs) _ I '^"^^^rjb'^s)'' ^iTj§j)) 

a-tn{s) J d'^h Gi„(6, s) 



where ip is the elementary multiplicity distribution and 
r{b,s) —< n >{b,s) / < n >{s) the elementary multiplic- 
ity function. 

In Ref. [2^, in the context of the geometrical picture, 
the elementary contact interaction process was based on 
e^e~ scattering data. In this approach we express 



r{s,b)^as)x]{s,b), 



where 



as) 



Jd%GUs,b) 

Jd%G,n{s,b)x]{s,b) 



and the power 7 is determined by fitting the average mul- 
tiphcity from e+e~ scattering data through a power law 
parametrization: 



V. EIKONAL MODELS 

It is expected that the eikonal function in the mo- 
mentum transfer space, xi^il): may be connected with 
some microscopic aspects of the underlying field theory 
(elementary interactions, form factors, structure func- 
tions) and, as mentioned before, it corresponds to a uni- 
tarized scheme connected with the experimental data. 
Eikonal models are characterized by different choices for 
x(s,g). In what follows we discuss our results and re- 
searches through two eikonal models (geometrical and 
QCD-based) 



A. Geometrical Model - Inelastic Channel 

In this subsection we review the description of pp mul- 
tiplicities distributions (inelastic channel) from models 
for the elastic channel in the context of the geometrical 
picture (contact interactions). 



< n >e+e-= ^[Vs]''. 

The elementary distribution ip{Z/r(b, s)) is repre- 
sented by a Gamma distribution and determined also by 
fits to e+e~ data. 

With inputs for Gi„(s,6) and/or Xi{s,b), obtained 
from fits to elastic scattering data, we have no free pa- 
rameter and the hadronic multiplicity distribution as 
function of Z and s may be inferred. 
- Elastic- channel inputs and results 

In Ref. (2^ we made use of three inputs from the 
elastic sector. Two are based on the Multiple Diffrac- 
tion Formalism, in which the eikonal in the momentum 
transfer space is expressed by 

X{b, s)^C j qdq.Mqb)GAGBf, (39) 

where Ga and Gb are the hadronic form factors, / the 
elementary (constituent - constituent) amplitude and C 
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does not depend on the transferred momentum. In this 
case we made use of the parametrizations used by Chou 
and Yang and also by Menon and Pimentel. Both present 
good descriptions of the experimental data in the elastic 
channel. The other input corresponds to the Short Range 
Expansion of the inelastic overlap function, introduced 
by Henzi and Valin, 



G,„(6,s) - P(s)exp{-&V4S(s)}fc(x,s), 

with k being expanded in terms of a short-range variable 
x = b exp{-(e6)V4B(s)}, i.e. 



N 

k{x, s) ^^6' 

n=0 



2n 



e cxp{l/2} 



With particular parametrizations excellent agreement 
with experimental data on pp and pp elastic scattering 
is achieved, allowing to infer the black-edge-large (BEL) 
behavior. 

In Ref. 01 a detailed discussion is presented on several 
variants from the elastic channel and parametrizations 
from e^e+ scattering. In particular, the results for the 
multiplicities distributions, with the BEL inelastic over- 
lap function, at -y/s = 52.6 and 546 GeV are displayed 
in Fig. 14, together with the experimental data. The 
prediction shows that the violation of the KNO scaling 
is well described. 



B. QCD-inspired models 

In this section we outline some research in course with 
the eikonal approach in connection with some QCD con- 
cepts. As we shall discuss the main point concerns the 
gluon-gluon contributions to the hadronic cross sections 
which we have investigated either from a dynamical gluon 
mass approach or by introducting the momentum scale 
in the gluonic distribution functions. After a review on 
the basic formalism we outline some aspects of both ap- 
proaches. 
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FIG. 14: Scaled multiplicity distribution for inelastic pp data 
at ^/s = 52.6 GeV and pp data at 546 GeV, compared with 
the model predictions {2^. 



1. Basic formalism 

The formalism was introduced by Afek, Leroy, Mar- 
golis and Valin '3^ and developed by several authors, 
including (for our purposes) Durand, Pi |3lj |. Block, Gre- 
gores, Halzen and Pancheri js j |. 

Originally, the point was to separate contributions 
from soft (S) and semi-hard (SH) inelastic processes by 
expressing 



where Ps is the probability of NO soft inelastic process 
and PsH the probability of NO semi-hard inelastic pro- 
cess. Therefore, that indicated an additive contribution 
in the eikonal: x{s,b) = xs(s,6) + Xsh(sJ>). 

In the recent version by Block et al. [S^l different ele- 
mentary contributions from quarks and gluons have been 
introduced: the gluon-gluon contribution comes from the 
parton model, the quark-quark from regge parametriza- 
tion and the quark-gluon by phenomenological inputs. In 
what follows we shortly review the main formulas. 

The normalization for the eikonal reads 



Gi„ei(s,6) 



PSPSH 

g-2_Re xs(s,b)g-2fle XSH{s,b) 



F{s,q) = ik / bdbJoiqb) 
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and the quark-gluon term as 



X{s,b) = Rex(s,6) +ilmx(s,6). 

For pp and scattering the crossing even and odd 
contributions are expressed by Xpp = +X~ a-^id Xpp = 
— x~ ■ The odd eikonal is assumed not to contribute 
at the asymptotic energies and is parametrized by 



Analyticity (generation of real and imaginary parts) 
for the even part is assumed as given by the prescription 



Re x^(s, b) + ilm x^(s, b). 



The even eikonal is expressed as a sum of three con- 
tributions, from quark-quark {qq), quark-gluon (qg) and 
gluon-gluon (gg) interactions, 



(Tqgis) = Cgg lOg — . 

So 

The gluon-gluon contribution is considered as the re- 
sponsible for the increase of the total cross section at 
the highest energies and is calculated through the parton 
model approach. 



, / dTFgg{T)agg{s), 

Jo 



(40) 



with 



^gg = / dxidx2fg{xi)fg{x2)S{T - X1X2), (41) 
Jo Jo 

where fg{xi) is the gluon distribution function, t = sis, 
and the symbol " denotes the elementary process. In |32l | 
the elementary cross section is given by 



X+{s, b) = Xqq{s, b) + Xqg{s, b) + Xgg{s, 5), 
which individually factorize in s and b, 

Xi]{s,b) = ia^j{s)w{b,^l.^■i), 

where i, j = q, g. 

The impact parameter distribution function for each 
process comes from convolution involving dipole form fac- 
tors (Chou-Yang Model): 

Wu{b,fiu) = f d^b'p,{\b'\p,{\b-b'\), 



p{b) =< G{q) > = < 



(l + g2/M')' 



where the angular brackets denote the symmetrical two- 
dimensional Fourier transform. Therefore, 



wuib) = lJi^[^^ubfK^i^lub), 

O iZTT 



and for i 7^ j it is assumed that 



The elementary cross sections for each process are in- 
troduced as follows. The quark-quark contribution is 
parametrized as a constant plus a Regge (even) term, 



^ggis) 



e{s - ml), 



(42) 



implying a cutoff mo for the particle production threshold 
and it is assumed the following simple parametrization 
for the gluon distribution function 



fgi^) 



N. 



{1-xf 



l (6-e)(5-6)...(l-e) 
2 5! 



(43) 



The model has 6 fixed parameters, mg, e, p-qq, fJ^gg, 
^■odd, Ols and 6 free parameters, determined from fits to pp 
and pp forward scattering data, namely atot{s), p{s) and 
B{s) above 15 GeV ^|. We have shown that the model 
applies only to forward and small momentum transfer 
regions .33]. 

In the next two sections we shall discuss two researches 
in course concerning the determination of the contribu- 
tion from gluon-gluon interactions, Eqs. (39-42). 



2. Dynamical gluon mass 

The possibility that the gluon propagator may be reg- 
ularized by a dynamically generated gluon mass [34j has 
recently provided important phenomenological descrip- 
tion of several processes |3^. The approach allows to 
calculate the contribution for the elementary gg cross 
section Eq. (42) and the main point is the association of 
the mass scale with the dynamical gluon mass. The basic 
ingredients are the expressions for the dynamical gluon 
mass, 
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FIG. 15; Description of the total cross sections through the 
QCD-based model with dynamical gluon mass mo = 500, 600 
and 700 MeV Isl. 



0-33(5) 



d\t\ 



with 



dan 



d\t\ 



■is J) 



97ra? , 3 



t 
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and by considering = \t\. 

The novel input concerns the updated determina- 
tions of the gluon distribution functions (CETEQ6), 
parametrized by means of Chebyshev polynomials. 
Presently, the implementation in the QCD-inspired ap- 
proach is being developed. 



VI. NONPERTURBATIVE QCD 



ln[(s + 4mg2)/A 



qcdJ 



-12/11 



and the associated running coupling constant 



as{s) 



where Po = 11 



An 



f3o\n {S + 4M^{s))/Al^^ 



3"/ 



and nf is the number of flavors. 



Preliminary tests with these contributions, in the context 
of the model described in the last subsection, have shown 
that the experimental data on CTtot(s), p{s) and B{s) are 
well described (13, as exemplified in Fig. (15). We are 
presently investigating the contributions from the other 
elementary processes, qq and qg. 



3. Momentum Scale 

Presently, we are attempting to improve the descrip- 
tions of the QCD-inspired models by taking into account 
the momentum transfer scale in the gluon distribution 
functions. The point is to replace the simple choice in 
Eq. (42), by distribution functions with the depen- 
dence, namely 



fgiXi) > fg{Xi,Q ). 

That can be implemented, following the approach by 
Durand and Pi jsj, by introducing the differential cross 
section 



As commented before the difhculties associated with 
high-energy soft processes arise from the fact that per- 
turbative QCD can not be applied and presently we 
do not know how to calculate even the elastic hadron- 
hadron scattering amplitudes from a pure nonperturba- 
tive QCD formalism. However, progresses have been 
achieved through the approach introduced by Landshoff 
and Nachtmann [s^l, developed by Nachtmann 38] and 
connected with the Stochastic Vacuum Model (SVM) (in- 
troductory reviews may be found in |39|). In particular, 
through this formalism and in some restricted kinematic 
conditions, it is possible to connect the gluon two-point 
correlation function with elementary (quark-quark) scat- 
tering amplitude. 

In this section we review the results we have obtained 
for these amplitudes with correlators determined from 
lattice QCD and also in the context of the Constrained 
Instantons. 



A. Stochastic Vacuum Model 

The approach has its origins in the attempts by Land- 
shoff and Nachtmann to connect soft high-energy pro- 
cesses with nonperturbative properties of the QCD vac- 
uum, as for example, the gluon condensate introduced by 
Shifman, Vainstein and Zakharov . In the first version 
jSJ] quarks couple with Abelian gluons. The non-Abelian 
version was developed by Nachtmann in the context of 
QCD and using the eikonal method for high energy in- 
teractions 38] . The scattering amplitude is calculated 
by means of a functional integral approach and is con- 
nected to a correlation function of two lightlike Wegner- 
Wilson loops. These correlation functions can be evalu- 
ated through the Stochastic Vacuum Model, in which the 
low frequency contributions to the functional integral of 
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QCD are described in terms of a stochastic process by 
means of a cluster expansion |4]|. The model incorpo- 
rates the gluon condensate concept and assumes that the 
correlation of two field strengths decreases rapidly with 
distance; due to an effective chromomagnetic monopole 
condensate, the QCD vacuum acts as a dual supercon- 
ductor. 

In this formalism the low frequencies contributions in 
the functional integral of QCD are described in terms 
of a stochastic process, by means of a cluster expansion. 
The most general form of the lowest cluster is t he g auge 
invariant two-point field strength correlator 0, |42| 

< FC^(a;)FD (y) 

+ d^{za5^p - Zp&^„)]{l - k)Di{z'^ /a)}, 

where z = x — y \s the two-point distance, a is a char- 
acteristic correlation length, k a constant, < FF > 
the gluon condensate and Nc the number of colours 
(C,D = 1, ...,N^ — 1). The two scalar functions D and 
Di describe the correlations and they play a central role 
in the application of the SVM to high energy scattering. 
Once one has information about D and Di, the SVM 
leads to the determination of the elementary quark-quark 
scattering amplitude, which constitutes important input 
for models aimed to construct hadronic amplitudes. The 
main formulas are as follows. 

The elementary amplitude / in the momentum transfer 
space is expressed in terms of the elementary profile 7 by 



After a two-dimensional integration, e{b) can be ex- 
pressed in terms of the correlation functions by 

e{b)^ei{b) + eii{b), (46) 

where 

/oo 
db'{b'-b)T^^[D{e)]{b'), (47) 
b 

euib) = (1 - n){g^FF)T^\-^D,{e)]{b). (48) 

YoyV = D oyV ^ Di we have V{k'^) = Ti[V{z'^)], where 
Tn denotes a n-dimensional Fourier transform. 

With the above formalism, once one has inputs for the 
correlation functions D{z) and Di{z), the elementary 
amplitude in the momentum transfer space, Eq. (43), 
may, in principle, be evaluated through Eqs. (44-47). It 
is important to stress that, as constructed, the formal- 
ism is intended for small momentum transfer [q^ ^ C'(l) 
GcV^) and asymptotic energies s ^ 00. Despite of these 
limitations, the investigation of soft high energy scatter- 
ing at the e nergies p resently available has led to satisfac- 
tory results lialia li ^ . 

B. Elementary Amplitudes 

In this section we review the results we have obtained 
from inputs for the above correlators from lattice QCD 
[4^ |4^ . We also comment the research in course in the 
semi-classical context of Instantons lig • 



/(q2) = / bdbJoiqbhib). (44) 
Jo 

In the Nachtmann approach the no-colour exchange 
parton-parton (loop-loop) amplitude can be written as 

where () means the functional integration over the gluon 
fields (the integrations are over the respective loop ar- 
eas), and w is a common reference point from which the 
integrations are performed. In the SVM by taking the 
Wilson loops on the light-cone the leading order contri- 
bution to the amplitude is given by 



1. Correlators from Lattice QCD 

Numerical determinations of the above correlation 
functions, in limited interval of physical distances, ex- 
ist from lattice QCD in both quenched approximation 
(absence of fermions) and full QCD (dynamical fermions 
included) 

With the procedure described above (see ^'^^ 
the calculational details), the elementary scattering am- 
plitude in the momentum transfer space can be deter- 
mined in numerical form. In order to obtain analyti- 
cal expressions, suitable for investigating distinct con- 
tributions and also for phenomenological uses, we have 
parametrized these numerical points through a sum of 
exponentials in q^: 



l{b)^r,e^b), (45) /(^ 

where rj is a constant depending on normalizations and 



The results are displayed in Fig. (16) from both quenched 
2 f J rr / IP f \ip f \\ approximation and full QCD, together with the corre- 

sponding exponential components. 



/ I dap,^dupaTr{Fp,^{x;w)Fpa(y;w)). 
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FIG. 16: Elementary amplitudes from quenched and lull 
QCD and the exponential components through parametriza- 
tion (48) [13. 



Our main conclusions are the following ji^: (1) the 
amplitudes decrease smoothly as the momentum transfer 
increases and they do not present zeros; (2) the decreas- 
ing is faster when going from quenched approximation to 
full QCD (with decreasing quark masses), and this effect 
is associated with the increase of the correlation lengths; 
(3) the dynamical fcrmions generate two contributions 
in the region of small momentum transfer, which are of 
the same order at ~ 1 GeV^ (only one contribution is 
present in the case of quenched approximation) . 

We understand that result (3) may suggest some kind 
of change in the dynamics at the elementary level, near 

^ 1 GeV^ and at asymptotic energies. If that is true, 
some signal could be expected at the hadronic level. One 
possibility is that this effect can be associated with the 
position of the dip (or the beginning of the "shoulder" ) in 



the hadronic (elastic) differential cross section data. The 
asymptotic condition embodied in our result indicates 
that ^ 1 GcV^ seems to be in agreement with the 
limit of the shrinkage of the diffraction peak, empirically 
verified when the energy increases in the region 23 GeV 
< ^ < 1.8 TeV. 

2. Correlators from the Instanton Approach 

By means of the stochastic vacuum formalism we also 
presently investigate the elementary amplitudes using 
correlators determined in the context of the constraint 
instanton approach, developed by A. Dorokhov and col- 
laborators • The basic picture is that of an instanton 
field dominating at small distances and decreasing expo- 
nentially at large distances in the physical vacuum. 

In Ref. [4^ , we make use of suitable parametrizations 
for the correlators and investigate the effects of the con- 
tributions from the short and long range correlations in 
the determination of the full correlator. Denoting those 
contributions as Di[z) and Dl{z), respectively, we in- 
troduce a dimensionless parameter a = rjgPc in terms of 
the driven parameter rjg and the size parameter pc '48'| . 
Since rjg is correlated with the relative contribution of 
each kind of correlator, we consider two extreme cases: 
1) equal contributions (weights 0.5 and 0.5), correspond- 
ing to a = 1.0; 2) almost pure instanton contribution 
(weights 0.99 and 0.01), corresponding to a = 0.1. 

In the lack of information on the long range compo- 
nent, and for our purpose, we consider parametrization 
in a Gaussian form j4q 

i?LW =exp{-(2/2.5)V}. 

For the short range case, a — 0.1, we introduce the 
parametrization, 

Di{z) = 0.7119cxp{-(2.403|2|)^} 
+ 0.2899exp{-(1.485|z|)2} 
- 9.456 X 10^3 exp{-1.277|z|}. 

The full correlator is then determined by 

D{z) = 0.99Di{z) + 0MDLiz). (50) 

For the long range case, a = 1.0, the parametrization 
takes the form 

Di{z) = 0.80084 exp{-(2.3025|z|)2} 

- 3.3846 X 10~2exp{-(0.97119|z|)2} 
+ 0.24225 exp{-2.7706|z|}, 

and for the full correlator we have 
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FIG. 17: Elementary amplitudes from the constrained instan- 
ton approach. 



D{z) = 0.5Dj{z) + 0.5Dl{z). (51) 

With the Eqs. (49) and (50) we can calculate the el- 
ementary amplitude through the steps indicated in Sec. 
VI. A. The results are displayed in Fig. (17). 

A central result is that if the contribution from the 
long range correlator is small, that is, an almost purely 
instanton case, the corresponding elementary amplitude 
presents a minimum. In terms of the associated differen- 
tial cross section, this implies a diffractive pattern in the 
momentum transfer space, a result already indicated in 
some phenomenological approaches ;45t . <46i | . In the case 
of equal weights the amplitude decreases monotonically 
with the momentum transfer. 

We conclude that in the context of the instanton ap- 
proach, the balance between the contributions of the 
short and long range correlators is a crucial point for the 
determination of the behavior of the elementary ampli- 
tudes. Further investigation along this line can bring new 
important insights on the connection between instanton 
correlators and the physical quantities which characterize 
the high-energy hadronic scattering. 



VII. PERSPECTIVES AND OUTLOOK 

In this section we outline some perspectives in the area 
of high-energy elastic hadron scattering from both exper- 
imental and theoretical point of views. Certainly, some 
ideas may be biased by our own knowledge and our own 
personal view. 



A. Experiments 

From the experimental point of view the perspectives 
are very optimistic due to the new generation of experi- 
ments with both accelerators and cosmic-ray observato- 
ries. Let us quickly summarize some projects in develop- 
ment. 

The upgrade of the Fermilab Tevatron machine, to- 
gether with upgrades and new devices in the CDF (Col- 
lider Detector Facility) and DO detectors, are going to 
provide improved investigations on pp collisions at ^/s 
2 TeV. Although the main purpose of the experiment 
concerns hard diffraction, it will be possible to investi- 
gate elastic scattering in both high and low regions, 
the slope parameter, total cross section and single diffrac- 
tion. Of topical importance for soft physics, the new de- 
termination of the total cross section shall possibly bring 
a solution for the puzzle represented by the discrepant 
results around 2 TeV (E710/E811 and CDF). 

At the Relativistic Heavy Ion Collider (RHIC) in 
the Brookhaven National Laboratory, pp collisions are 
presently being investigated at energies never reached be- 
fore: Vs : 50 - 500 GeV. The experiment "Total and 
Differential Cross Sections and Polarization Effects in pp 
Elastic Scattering at RHIC" (pp2pp) plans to investigate 
both elastic scattering and diffraction dissociation (single 
and double), in addition to spin effects. This will provide 
the first opportunity for direct comparison between pp 
and pp scattering at the highest collider energies. 

Although in a bit longer term, at the CERN Large 
Hadron Collider (LHC), the TOTEM experiment (Total 
Cross Section, Elastic Scattering and Diffraction Disso- 
ciation at the LHC) is specifically planned to study soft 
diffractive physics in pp collisions at ~ 16 TeV. In 
particular, diffraction dissociation, total cross section and 
elastic scattering at large values of the momentum trans- 
fer will be investigated, up to = 10 — 15 GeV^. That 
will certainly allow discrimination and selection of vari- 
ous models and approaches, giving fundamental informa- 
tion at large momentum transfer; for example, showing 
the existence or not of structures. Moreover, this experi- 
ment will probably provide a final answer on the possible 
differences between pp and pp total cross sections and the 
correct power 7 in the In'' s dependence of atot[s). 

The most energetic event detected in cosmic-ray ex- 
periments had Eiab = 3 X 10^" eV, corresponding to an 
energy of 50 Joules! Goals of the Auger project are the 
measurement of arrival direction and the energy and mass 
composition of cosmic rays above Eiab = 10^^ eV. For 
pp collisions this means ^/s above 140 TeV, nearly 10 
times the LHC energy. In addition to the astrophysical 
importance of the experiment, the measurement of the 
longitudinal development of showers will provide severe 
tests on hadronic interaction models. As a consequence, 
among others, the puzzles concerning the extraction of 
pp cross section from p — air production cross section 
may receive better insights, allowing more precise deter- 
minations and at energies possibly never to be reached 
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by accelerator machines. 

B. Theory 

Elastic hadron scattering (and soft difFractive processes 
in general) is a long distance phenomena and there- 
fore we expect and look for a theoretical treatment via 
non-perturbative QCD. Despite all the difficulties men- 
tioned along this manuscript, we understand that two 
approaches deserve special attention. 

One of them concerns the approach by Nachtmann and 
the Stochastic Vacuum Model (Sec. VI. A) jHlllil- Al- 
though under restrictive kinematic conditions (momen- 
tum transfer of the order or below 1 GeV^ , and asymp- 
totic energies, s —^ oo) the formalism has provided inter- 
esting results in the investigation of the physical quan- 
tities that characterize the elastic pp and pp scattering 
[4^ li^ . in special the works by Ferreira and Pereira, 
connecting experimental observables and QCD parame- 
ters "43] . Attempts to implement the energy dependence 
in pure QCD grounds may be an important task for the 
near future. 

The other approach is associated with evidences for 
finite gluon propagator and running coupling in the in- 
frared region and that is the case in some classes of solu- 
tions of the nonperturbative Schwinger-Dyson equations. 
In particular, in the solution proposed by Cornwall [33 | 
the gluon acquires a dynamical mass leading to a freez- 
ing of the coupling constant in the infrared region. As 
referred before, Natale and collaborators have dis- 
cussed several phenomenological tests, reaching interest- 
ing results which have permitted the development of the 
formalism and the selection of adequate basic inputs. 
That opens a new way to investigate long distance phe- 
nomena with a finite calculational approach. 

We also understand that the connections between soft 
and semi-hard processes, typical of QCD-based models in 
the eikonal context, may bring new insights for the de- 
velopment of adequate calculational schemes in the non- 
perturbative treatment of high-energy elastic collisions. 

VIII. SUMMARY AND FINAL REMARKS 

Despite its simplicity, elastic hadron scattering consti- 
tutes a topical problem in high-energy Physics. Although 
the bulk of experimental data can be efficiently described 
in different phenomenological contexts, we are still facing 
the lack of a treatment and of a reasonable understanding 
of these processes based exclusively on QCD. 

Our main strategy in investigating elastic scattering 
has been to look for descriptions based on the high-energy 
principles and theorems from Axiomatic Quantum Field 
Theory and, simultaneously, attempting to extract "em- 
pirical" information from all the experimental data avail- 
able. Tests of discrepant data and their influence on the 
extracted information play a central role. In that way 



we hope to get feedbacks for theoretical development in 
nonperturbative and semi-hard QCD. We can summarize 
our main recent results as follows. 

In the context of the Analytic Approach, we have in- 
vestigated the effects of discrepant experimental informa- 
tion on the total cross sections in both accelerator and 
cosmic-ray energy regions. By means of analytical fits, we 
have obtained extrema bounds for the soft Pomeron in- 
tercept, namely ^^^''''(O) = 1.109 and a5g'"'=''(0) = 1.081. 
We have also obtained novel constrained bounds for the 
intercept from spectroscopy data (fitted Regge trajecto- 
ries from Chew-Fautschi plots) and extended the analysis 
to several reactions. That information on the Pomeron 
bounds may be important for phenomenological develop- 
ments and projects for new experiments. We have also 
shown that the presence of the Odderon in the real part 
of the elastic hadronic amplitude is not forbidden by the 
bulk of experimental data on atot and p. In particu- 
lar, the fit with the Kang-Nicolescu parametrization has 
indicated a crossing in p{s), with ppp becoming greater 
than Ppp above y/s = 70 GeV. That parametrization pre- 
dicts ppp{y/s = 200GeV) = 0.134±0.005 (RHIC regions). 
Detailed investigation on the applicability of DDR have 
shown that, once the subtraction constant is used as a 
free fit parameter, the DDR is equivalent to the IDR 
with finite lower limit (sq = 2m^). That result was ob- 
tained for the class of entire functions in the logarithm 
of the energy (typical of analytic models) . 

In the context of Model Independent Analyses, we have 
investigated the correlations between the experimental 
data on total cross section and the slope parameter. 
The parametrization introduced is based on the empir- 
ical behavior of these quantities and extrapolations to 
cosmic-ray energies may be useful in the determination of 
proton-proton total cross sections from proton-air cross 
sections. By means of a novel model independent fit 
procedure to the differential cross section data, we have 
found statistical evidence for eikonal zero in the momen- 
tum transfer space and that the position of the zeros 
decreases as the energy increases. The zero position 
shows agreement with the result recently obtained for 
the electromagnetic form factor and inferred from elas- 
tic electromagnetic e~p scattering (polarization transfer 
experiments). Since our analysis concerns only hadronic 
interactions the results may bring new insights in the in- 
vestigation of electromagnetic and hadronic form factors. 
We are also treating analytic fits with free parameters 
depending on the energy so as to develop a model inde- 
pendent predictive approach. 

In the context of Eikonal Models we have obtained 
connections between the elastic and inelastic channel 
by means of the geometrical picture, and correlating pp 
and pp scattering with contact interactions simulated by 
e^e+ distributions and multiplicities. With the class of 
QCD-based or QCD-inspired models we have developed 
novel gluonic contributions by means of two approaches. 
One is based in solutions of the Schwinger-Dyson equa- 
tions characterized by the dynamical gluon mass. The 
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other one is intended to take into account the Q^- scale 
in updated gluonic distribution functions. Certainly, the 
two approaches are not independent, and we presently 
investigate their simultaneous implementation. 

In the context of the Stochastic Vacuum Model, we 
have obtained novel results for the elementary scatter- 
ing amplitudes, making use of correlators determined ei- 
ther from lattice QCD or from constrained instantons. 
In both cases the elementary amplitudes present no ze- 
ros (change of sign in the momentum transfer space). 
In the context of eikonal models, in which the eikonal 
is expressed in terms of the elementary amplitudes and 
form factors, Eq. (38), this result corroborates the in- 
terpretation of the zero in the hadronic form factor, and 
therefore, its dependence on the energy. 



Acknowledgments 



I am thankfuU to Y. Hama and F.S. Navarra for the 
invitation to present this review. For fruitful comments 
and discussions along these Workshops I am particularly 
grateful to E. Ferreira, Y. Hama and T. Kodama. I am 
deeply thankful to those who have developed research 
programs under my supervision and are co-authors of all 
the works reviewed here: R.F. Avila, P.C. Beggio, S.D. 
Campos, P.A.S. Carvalho, E.G.S. Luna, A.F. Martini, 
J. Montanha, and D.S. Thober. I am also thankful to 
A. Di Giacomo, P. Valin, A. Dorokhov, A. A. Natale and 
R.C. Rigitano for valuable discussions. This work was 
supported by FAPESP (Contract N. 00/04422-7). 



[1] Y. Dokshitzer, Phil. Trans. Soc. London A395, 309 

(2001) . 

[2] M. M. Block and R. N. Cahn, Rev. Mod. Phys. 57, 563 
(1985). 

[3] G. Matthiae, Rep. Prog. Phys. 57, 743 (1994). 

[4] E. Engel, Nucl. Phys. B (Proc. Suppl.) 82, 221 (2000). 

[5] E. Predazzi, in International Worksfiop on Hadron 
Physics 98, edited by E. Ferreira, F.F.S. Cruz and S.S. 
Avancini (World Scientific, Singapore, 1999) p. 80 - 175; 
V. Barone and E. Predazzi, High-Energy Particle Diffrac- 
tion (Springer- Verlag, Berlim, 2002). 

[6] S. Donnachie, G. Dosch, P. Landshofi', O. Nachtmann, 
Pomeron Physics in QCD, Cambridge University Press, 
Cambridge, 2002. 

[7] M.J. Menon, "Introduction to Soft Diffraction: Some Re- 
sults and Open Problems", LISHEP 2002, Session B: 
Advanced School on High Energy Physics (2002). 

[8] A.F. Martini and E. Predazzi, Phys. Rev. D 66, 034029 

(2002) . 

[9] K. Hagiwara et al, Phys. Rev. D 66, 010001 (2002), Re- 
view of Particle Physics, Particle Data Group. Full data 
sets are available at http://pdg.lbl.gov 

[10] R. F. Avila, E. G. S. Luna, and M. J. Menon, Phys. Rev. 
D 67, 054020 (2003). 

[11] R.J. Eden, Rev. Mod. Phys. 43, 15 (1971); S.M. Roy, 
Phys. Rep. 5, 125 (1972); J. Fischer, Phys. Rep. 76, 157 
(1981); P. Valin, Phys. Rep. 203, 233 (1991). 

[12] M. J. Menon, A. E. Motter, and B. M. Pimentel, Phys. 
Lett. B 451 , 207 (1999). 

[13] R. F. Avila, E. G. S. Luna, and M. J. Menon, Braz. J. 
Phys. 31, 567 (2001). 

[14] R.F. Avila, E.G.S. Luna, and M.J. Menon, in Proceedings 
of the Vni International Workshop on Hadron Physics 
2002, edited by C.A.Z. Vasconcellos et al. (World Scien- 
tific, Singapore, 2003) p. 335. 

[15] E.G.S. Luna and M.J. Menon, Phys. Lett. B 565, 123 

(2003) . 

[16] R.F. Avila and M.J. Menon, Nucl. Phys. A, 744, 249 

(2004) . 

[17] A.F. Martini, M.J. Menon, Phys. Rev. D 56, 4338 (1997). 

[18] E.G.S. Luna, M.J. Menon, in New States of Matter m 
Hadronic Interactions, PASI, AIP Conf. Proc. Vol. 631, 
edited by H.-Thomas Elze et al. (AIP. New York, 2002) 



pp. 721-722. 

[19] A. Donnachie and P.V. Landshoff, Z. Phys. C 2, 55 
(1979); Phys. Lett. B 296, 227 (1992); Phys. Lett. B 387, 
637 (1996). 

[20] K. Kang and B. Nicolescu, Phys. Rev. D 11, 2461 (1975); 
L. Lukaszuk and B. Nicolescu, Lett. Nuovo Cimento 8, 
405 (1973). 

[21] E.G.S. Luna, M.J. Menon, J. Montanha, Braz. J. Phys. 
34, 268 (2004). 

[22] E.G.S. Luna, M.J. Menon, J. Montanha, Nucl. Phys. A 

745, 104 (2004). 
[23] P.A.S. Carvalho, A.F. Martini, M.J. Menon, Eur. Phys. 

J. C 39, 359 (2005). 
[24] P.A.S. Carvalho and M.J. Menon, Phys. Rev. D 56, 7321 

(1997). 

[25] A.F. Martini, M.J. Menon, J. Montanha, "Zeros in the 
Electromagnetic and Hadronic Form Factors", in IX 
Hadron Physics and VIII Relativistic Aspects of Nuclear 
Physics - A Joint Meeting on QCD and QGP, edited by 
M.E. Bracco, M. Chiapparini, E. Farreira, T. Kodama, 
AIP Conf. Proc. Vol. 739 (AIP, New York, 2004) pp. 578 

- 580. 

[26] P.A.S. Carvalho, A.F. Martini, M.J. Menon, A.F. Mot- 
ter, in International Workshop on Hadron Physics 98, 
Florianopolis, 1998, edited by E. Ferreira, F.F. de Souza 
Cruz and S.S. Avancini (World Scientific, Singapore, 

1999) p. 326. 

[27] R.F. Avila, S.D. Campos, M.J. Menon, and J. Montanha, 
"Analytical Fits to Hadron-Hadron Differential Cross 
Section Data at the Diffraction Peak", in IX Hadron 
Physics and VIII Relativistic Aspects of Nuclear Physics 

- A Joint Meeting on QCD and QGP, edited by M.E. 
Bracco, M. Chiapparini, E. Farreira, T. Kodama, AIP 
Conf. Proc. Vol. 739 (AIP, New York, 2004) pp. 529 - 
531. 

[28] A.F. Martini, M.J. Menon, J. Montanha, Braz. J. Phys. 
34, 263 (2004). 

[29] P.C. Beggio, M.J. Menon, and P. Valin, Phys. Rev. D 
61, 034015 (2000); in International Conference on Elastic 
and Diffractive Scattering - Blois VIII, edited by V.A. 
Petrov and A.V. Prokudin (World Scientific, Singapore, 

2000) pp. 371 - 376. 

[30] Y. Afek, C. Leroy, B. Margolis, P. Valin, Phys. Rev. Lett. 



24 



45, 85 (1980). 

[31] L. Durand and H. Pi, Phys. Rev. Lett. 58, 303 (1987); 
Phys. Rev. D 38, 78 (1988); Phys. Rev. D 40, 1436 (1989). 

[32] M.M. Block, E.M. Gregores, F. Halzen, G. Pancheri, 
Phys. Rev. D 60, 054024 (1999). 

[33] E.G.S. Luna, A.F. Martini, and M.J. Menon, in Pro- 
ceedings of the VIII International Workshop on Hadron 
Physics 2002, edited by C.A.Z. Vasconcellos et aL (World 
Scientific, Singapore, 2003) p. 295. 

[34] J. M. Cornwall, Phys. Rev. D26 (1982) 1453; J. M. Corn- 
wall and J. Papavassiliou, Phys. Rev. D40 (1989) 3474, 
D44 (1991) 1285. 

[35] M. B. Gay Ducati, F. Halzen and A. A. Natale, Phys. 
Rev. D48, 2324 (1993); F. Halzen, G. Krein and A. 
A. Natale, Phys. Rev. D47, 295 (1993); A. C. Aguilar, 
A. Mihara and A. A. Natale, Phys. Rev. D65, 054011 
(2002), Int. J. Mod. Phys. A19, 249 (2004); A. C. 
Aguilar, P. S. R. da Silva and A. A. Natale, Phys. Rev. 
Lett. 90, 152001 (2003); 

[36] E.G.S. Luna, A.F. Martini, M.J. Menon, A. Mihara, A.A. 
Natale, "Preliminary results on the influence of a gluon 
mass in hadronic scattering" in IX Hadron Physics and 
VIII Relativistic Aspects of Nuclear Physics - A Joint 
Meeting on QCD and QGP, edited by M.E. Bracco, M. 
Chiapparini, E. Faxreira, T. Kodama, AIP Conf. Proc. 
Vol. 739 (AIP, New York, 2004) pp. 572 - 574. 

[37] P.V. Landshoff, O. Nachtmann, Z. Phys. C 35, 405 
(1987). 

[38] O. Nachtmann, Ann. Phys. 209, 436 (1991). 

[39] O. Nachtmann, "High-Energy Collisions and Nonpertur- 
bative QCD", hep-ph/96093; O. Nachtmann, in Inter- 
national Workshop on Hadron Physics 2000, edited by 
F.S. Navarra, M.R. Robilotta, G. Krein (World Scientific, 
2001) p. 152-216; E.G. Dosch. in Intenuitional Workshop 
on Hadron Physics 96, edited by E. Ferreira et al. (World 
Scientific, 1997) pp. 169-219. 



[40] M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. 

Phys. B 147, 385, 448, 519 (1979). 
[41] E.G. Dosch, Phys. Lett. B 190, 177 (1987); E.G. Dosch, 

Yu. A. Simonov, Phys. Lett. B 205, 339 (1988); Yu. A. 

Simonov, Nucl. Phys. B 307, 512 (1988); 204, 67 (1989); 

A. Kramer, E.G. Dosch, Phys. Lett. B 252, 669 (1990); 

272, 114 (1991). 
[42] E.G. Dosch, E. Ferreira and A. Kramer, Phys. Rev. D 

50, 1992 (1994); Phys. Lett. B 318, 197 (1993); 289, 

153 (1992). 

[43] E. Ferreira and F. Pereira, Phys. Rev. D 55, 130 (1997); 
Phys. Rev. D 56, 179 (1997); Phys. Lett. B 399, 177 

(1997); Phys. Rev. D 59, 014008 (1998); Phys. Rev. D 
61, 077507 (2000); E. Ferreira, Braz. J. Phys. 30, 282 
(2000). 

[44] E.R. Berger, O. Nachtmann, Eur. Phys. J. C 7, 459 
(1999). 

[45] M.J. Menon, Nucl. Phys. B (Proc. Suppl.) 25, 94 (1992); 
Phys. Rev. D 48, 2007 (1993); Can. J. Phys. 74, 594 

(1996). 

[46] A.F. Martini, M.J. Menon, D.S. Thober, Phys. Rev. D 

57, 3026 (1998). 
[47] A.F. Martini and M.J. Menon, Phys. Lett. B 570, 53 

(2003). 

[48] A. Dorokhov. A.F. Martini, M.J. Menon, in XX Brazil- 
ian National Meeting on Particles and Fields (1999), 
www.sbf 1 . if . usp . br/ eventos /enfpc/xx/procs/ res68/ . 

[49] A. Di Giacomo and E. Panagopoulos, Phys. Lett. B 
285 (1992) 133; A. Di Giacomo, E. Meggiolaro, and E. 
Panagopoulos, Nucl. Phys. B 483 (1997) 371; M. D'Elia, 
A. Di Giacomo, E. Meggiolaro, Phys. Lett. B 408 (1997) 
315. 

[50] A.E. Dorokhov, S.V. Esaibegian, A.E. Maximov, S.V. 
Mikhailov, Eur. Phys. J. C 13, 331 (2000). 



